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ABSTRACT 

We obtained X-ray XMM-Newton observations of the open cluster Berkeley 87 and the massive star- 
forming region (SFR) ON 2. In addition, archival infrared Spitzer Space Telescope observations were used 
to study the morphology of ON 2, to uncover young stellar objects, and to investigate their relationship 
with the X-ray sources. It is likely that the SFR ON 2 and Berkeley 87 are at the same distance, 1.23 
kpc, and hence are associated. The XMM-Newton observations detected X-rays from massive stars in 
Berkeley 87 as well as diffuse emission from the SFR ON 2. The two patches of diffuse X-ray emission are 
encompassed in the shell-like H n region GAL 75.84+0.40 in the northern part of ON 2 and in the ON 2S 
region in the southern part of ON 2. The diffuse emission from GAL 75.84+0.40 suffers an absorption 
column equivalent to Ay « 28 mag. Its spectrum can be fitted either with a thermal plasma model at 
r;30MK or by an absorbed power-law model with y«-2.6. The X-ray luminosity of GAL 75.84+0.40 is 
Lx~l x 10 32 erg s _1 . The diffuse emission from ON 2S is adjacent to the ultra-compact H n (UCH n) region 
Cygnus 2N, but does not coincide with it or with any other known UCH n region. It has a luminosity of 
Lx~6 x 10 31 erg s -1 . The spectrum can b e fitted with an absorbed power-law model with y«— 1.4. We 
adopt the view of iTurner & Forbes dl982l) that the SFR ON 2 is physically associated with the massive 



star cluster Berkeley 87 hosting the WO type star WR 142. We discuss different explanations for the 
apparently diffuse X-ray emission in these SFRs. These include synchrotron radiation, invoked by the 
co-existence of strongly shocked stellar winds and turbulent magnetic fields in the star-forming complex, 
cluster wind emission, or an unresolved population of point sources. 

Subject headings: H II regions - open clusters and associations: individual (Berkeley 87) - stars: early-type - stars: 
winds, outflows - X-rays: ISM - X-rays: stars 



1. Introduction 

Recent observational advances in X-ray astro- 
physics have led to a new high-energy perspective 
on the interstellar medium and star-forming regions 
(SFRs). The X-ray point sources in SFRs comprise 
pre-main sequence (PMS) stars as well as massive 
stars. Sometimes, X-ray emission from deeply em- 
bedded young stellar objects ( YSO) is observed, al beit 
such observations remain rare dPravdo et al. I l2009h . It 
is now firmly established that the shocked winds of 
young massive OB stars contri bute to the heating of 
interstellar matter up to JlOMK (iTownsley et al .120031 



Gudeletal ]l2008h . In older massive stellar aggregates, 



the combined action of stellar winds and supernova ex- 
plosions results in powerful cluster winds an d associ- 
ated superbubbles ( ChevalieL&£leggl 1985). Higher 
temperatures (up to 100MK) are observed in these ob- 
jects. Hot, X-ray-emitting gas is present around mas- 
sive stars and fills the large volumes of star clusters, 
and, in some cases beyond. 

There is a small group of SFRs where hard, some- 
times non-thermal X-ray emission is observed. In 
cases when strong fluorescent lines are seen in the 
spectra, the emission is explained by the presence of 
a recent supernova remnant (SNR) and its interaction 
with t he cool dense mat erial of a nearby molecular 
cloud (ITakagi et al.ll2002l) . However, in some promi- 
nent SFRs, such as RCW 38, the non-thermal X-ray 
emission cannot be easily explain ed unless the pres - 
ence of magnetic fields is assumed ( W olk et al 1 l2002h . 
Magnetic fields, along with tur bulence, play 



an im portant role in star formation (ICrutcher et al 



2009|) . Magnetic fields were directly measured in 



M 17 dBrogan & Trolandl200l1) and Orion dSchleuning 
>98l) . a nd are predict ed to be a common feature 
in SFRs dFerlandll2009l) . If accelerated particles are 



present in the same volume, synchrotron emission 
would naturally occur. 

A region within a few arcminutes of the OH maser 
ON 2 in Cygnus X is an established site of ongoing 
massive star formation. Following the literature we 



refer to this whole S FR as ON 2 dDent et alj 11988 



Shep herd et al.l 1 1997h . This region is located within 



the well-studied massive star cluster Berkeley 87. 

One of the most interesting members of Berkeley 87 
is the WO type star WR 142. Only three stars of this 
spectral type are known in the Galaxy, and the WO 
star in Berkeley 87 is the closest among them. Stars 
of this spectral type represent the evolutionary stage 



immediately before the supernova or y-ray burst ex- 
plosion, and drive the fastest stellar winds among all 
stars. Berkeley 87 traditionally attracts the attention of 
high-energy astrophysics as a potential site of particle 
acceleration. Therefore, SFR ON 2 provides an ideal 
laboratory to study the interactions between an active 
star-forming region and the massive star feedback. 

In this paper we present our XMM-Newton obser- 
vations of the field encompassing ON 2 and Berke- 
ley 87, and provide their analysis and interpretation. 
The paper is organized as follows. In Section [2] we 
introduce the SFR ON 2, and discuss its distance and 
relation to the star cluster Berkeley 87. The member- 
ship and evolution of Berkeley 87 is treated in Sect.[3] 
The analyses of XMM-Newton and Spitzer observa- 
tions are presented in Sect. [4] In Sect. [5] we briefly 
address the distribution of point sources detected by 
Spitzer. In Sect. [6] we consider X-rays from massive 
stars in Berkeley 87. Section|7]is devoted to the X-ray 
emission from the Hn region GAL 75.84+0.40. Sec- 
tion [8] concentrates on the observed properties of X- 
ray emission from ON 2S, the southern part of ON 2, 
while discussion on its origin is presented in Sect. [9] A 
comparison to other SFRs is presented in Sect.Qj)] and 
conclusions are drawn in Sect.QT| In the Appendix we 
briefly review the suggestions in the literature about 
the possible identification of y-ray sources in Berke- 
ley 87. 

2. The massive star-forming region ON 2 

FigureQ] shows the composite DSS optical and 
Spitzer IR image of SFR ON 2 and star cluster Berke- 
ley 87. A large molecular cloud complex, Onsala2C, 
observed in CO, is intimately associated with the 
whole ON 2 comple x dMatthews et al 1 19861). In a mor- 
phological model bv lTurner & Forbe sl dl982h based on 
studies of the reddening and obscuration, Berkeley 87 
sits on the western edge of the heavily obscured SFR 
ON 2 and the associated giant molecular cloud On- 
sala2C. Figure|2] gives a schematic representation of 
the morphology of the ON 2 region. The center of 
the massive star cluster Berkeley 87 is in the middle 
of the figure. The molecular cloud occupies the up- 
per left quadrant. The sites of active star formation, as 
highlighted by compact H n regions, are located at the 
edges of a molecular cloud. 

A Spitzer IRAC image of the r egion is shown in 
the right panel of Fig.[T] Following [Dent etalj dl988l) 
we will distinguish between northern and southern 
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Fig. 1. — Left panel: Combined optical and IR image of Berkeley 87 and ON 2. The POSS2/blue image (Copyright 
Second Palomar Sky Survey 1993-1995 by the California Institute of Technology) is shown in blue, the Spitzer IRAC 
channel 1 (3.6 pm) is in green, and channel 4 (8 pm) in red. T he large white circle represents the outer cluster boundary 
as determined from star counts by iTurner & Forbesl (119821) . The center of the circle is at 20 h 21 m 37 s , +37° 24' 37" 
(J2000), the radius is 8'. The Berkeley 87 cluster is fully within the XMM-Newton field-of-view (30')- The approximate 
extend of the massive star-forming region ON 2 is shown by the dashed square. The region withing the square is 
enlarged in the right panel. Right panel: Combined IR Spitzer IRAC (3.6 pm blue, 4.5 pm green, 8.0 pm red) image of 
the massive star-forming region ON2. The H n regions described in the text are identified by arrows. The approximate 
location of the molecular cloud Onsala2C is shown as a dashed circle. Image size is «12'xl 1'. North is up, and east 
is left. 



Hn regions, and introduce the notations ON2N and 
ON2S. The northern part, ON2N, comprises the 
Hn regions GAL 75.84+0.40 and GAL 75.84+0.36. 
ON2S contains the Hn regions Cygnus2N (alias 
G75.78+0.34, Onsala2N), Onsala2 (alias [HLB98] 
Onsalal30), [L89b] 75.767+00.344, and [HLB98] 
Onsalal29 (Fig.Q). The UCH ii regions Cygnus 2N 
and Onsala2 are separated by on ly » 2'.' 5 On- 



sala 2, Cygnus 2N, and Onsala 129 (Palagi et al.ll 1 993 
Hofner & Churchwelll[l996llShepherd et al.lll997h 



She pherd et al. I (1 19971) studied ON 2S in molecular 



The Simbad data-base gives aliases of Onsala 2 as 
[WAM82] 075.77+0.34, OH 75.8+0.3, [HLB98] Onsala 130, 
and [PCC93]414. The coordinates of Onsala 2 from Simbad are 
20 h 21 m 43 s 8, +37° 26' 39" (J2000). However, the coordinates of 
[WAM82] 075.77+ 0.34 are 20 h 21 m 4H31, +37° 25' 53.5" (J2000) 
Jwinketal.|[l982l) . T his object is »54" away from Onsala 2. The 
positional accuracy of IWink et al.U 1 9821 ) is typically 5". Therefore, 
[WAM82] 075.77+0.34 is not Onsala 2. 



lines and the 3 mm continuum. They detected three 
deeply embedded YSOs in Cygnus 2N, one of which 
is likely to be the driving engine of a molecular out- 
flow. From the dynamical timescales of the outflows 
(~30-50kyr) and the high luminosities inferred for the 
YSOs (5000, 430, and 330 L ) they suggest that a near- 
simultaneous massive star formation event occurred in 
this region ~ 10 4 yr ago. This su pport similar conclu- 
sions made bv lDent et alj (119881) . 



2.1. 



Distance to ON 2 and its relation to Berke- 
ley 87 



ON 2 is located in the Cygnus region, and thus 
we are looking tangentially to the Orion local spi- 



ral arm and observe numerous ob 



distances (lUvaniker et al.l |2Q01 ). 



ects at different 



Reifenstein et al 



(1 19701) estimate the dist ance to ON 2 as 5.5 kpc, while 
Turner & Forbesl (119821) argue that ON 2 is located at 
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675.84.0.40 




Fig. 2. — Sketch of th e possible morphology o f ON 2 
based on Fig. 6 from iTurner & Forbes! (1982). The 
black circle is the same as in Fig.[T] The shaded area 
follows contour lines joining regions of similar space 
reddening (units of 0.01 mag in Eb-v)- The pink el- 
liptical regions represent approximately the patches of 
hard X-ray emission. The blue dots represent stars in 
Berkeley 87, the red spots represent UCHn regions. 
Not to scale. 



the same distance as the cluster Berke ley 87 and is 



physic ally connected with this cluster. iMassev et al 



(12001b obtained new photometric measurements of 
stars in Berkeley 87 and der ived a distance of d ~ 



1600pc. ITurner et al.1 (120061) augmented their previ 



ous studies of Berkeley 87 by a larger number of stars, 
and derived d » 1230 + 40 pc, which we adopt here. 



Turner & Forbesl (119821) notice that a trunk-like 



zone of heavy obscuration seen in Fig. [2] corresponds 
spatially with a CO clouds bel onging to the Cygnus X 
complex. In a recent study ISchneider et al.l (2007) 
show that the UV radiation from the clusters within 
the Cygnus OBI association, including Berkeley 87, 
affects the molecular cloud complex in the Cygnus X 
south region, and note that a distance between 1 . 1 and 
1 .3 kpc is favored from O stars spectroscopy. This sup- 
ports the likely association between Berkeley 87 and 
ON 2. 

3. The open star cluster Berkeley 87 

Berkeley 87 is a relatively sparse, moderately red- 
dened cluster, whose members represent some of 



the key stages in the evolution of massive stars. 
The most intriguing is a WO-type star, WR 142. 
Analysis of its X-ray observations was presented in 



Oskinova et al 



The brightest member of 
Berkeley 87, HD 229059, lies slightly off the cluster 
core. It is a binary system with a B1.5Iap and a lower- 



lumin osi ty late-O or B0 star companion ( Negueruelal 



2004- iMathysl (119871) suggests that this star is a 



blue straggler. Close to the clust er center is t he pe 



culiar variable star V439Cyg. iNegueruelal (120041) 
identi fy this star as havi ng B1.5Ve spectral type, 
while IMassev et alJ (1200 II) suggest a B[e] classifica- 
tion. The bright red s upergiant BCCyg of M3.5Ia 
type (ITurner et al.l 120061) is located at the north of the 
cluster, close to the Hn region GAL 75.84+0.40. A 
possible s pectroscopic binary BD+36°4032 (Q8 .5 V or 
08.5111) dNegueruelall2004l: IMassev et al.ll200l"l) is lo- 
cated immediately south of ON 2S. 

The history of Berkeley 87 is not well understood. 
The ages of OB stars, the WO star, and the red super- 
giant star are in apparent disagreement with an evolu- 
tion within a nearly coeval cluster. 



Vanbeveren et al.l d 19981) found that 27% of the 



young Galactic massive clusters contain WR and red 
supergiant (RSG) members, and that such clusters 
must be older than 4Myr. Moreover, WC/WO stars 
can co-exist with RSG stars onl y during a short time 



interval of ~ a few x 10 5 yr. IMassev et alJ (1200 ll) 
studied Berkeley 87 among other Galactic open clus- 
ters. They suggested that OB stars in Berkeley 87 were 
formed within the time span of <1 Myr and are coeval 
with WR 142 at 3.2 Myr. Using evolutionary tracks of 
Schaller et al. (1992), they estimate the initial mass of 
WR 142 as Mj = 70 M Q . In the meantime, the in itial 
mass of BCCyg is S25M dLevesque et alJl2005h . A 
star of such mass must be older than 6. 4 Myr when it 
reaches RSG stage dSchalleret al.lll992b . 

This apparent contradiction can be resolved with 
stella r evolutionary models t hat account for rota- 
tion dMeynet & Maeden 120051) . We have recently 
found indications tha t WR 142 may be a fast rotator 
dOskinova et al.ll2009b . In fast rotating stars, the RSG 
stage will occur earlier, perhaps even during the H- 
burning stage, while the WO stage occurs later com- 
pared to the non-rotating models. Therefore, the si- 
multaneous presence of a RSG and a WO star can be 
explained if Berkeley 87 is ~ 4-6 Myr old. 

Assuming that the most massive star in Berkeley 87 
has an initial mass higher than 80 M and the univer- 
sal initial mass function (IMF), we estimate that about 
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Fig. 3.— Combined XMM-Newton MOS1, MOS2, and PN image of the Berkeley 87 in the 0.25 to 12.0keV band 
(red=0.25-1.0keV, green= 1. 0-2.5 keV, blue=2.5-12.0keV). Image size is *30'x30'. North is up, and east is left. 



30 stars with initial mass more than 1OM should be 
present in Berkeley 87. This estimate is consistent with 
the 22 know n massive stars in B erkeley 87 given the 
uncertainties ( Massev et alj|2001 ). In this case, the to- 
tal mass of Berkeley 87 should be « 1200 M e . The 
number of low mass stars with masses between 0.5 M 
and 3 M Q can be estimated as a 2000. 

4. XMM-Newton and Spitzer observations of 
Berkeley 87. 

Berkeley 87 was observed by XMM-Newton dur- 
ing two consecutive satellite orbits 0. The data were 
merged and analyzed using the software sas 8.0.0. Af- 
ter the high-background time intervals have been re- 
jected, the combined exposure time of all detectors 
was « lOOks. We followed the standard procedure 
for the source detection, setting the minimum detec- 



2 ObsId 0550220 10 1 , Obsld 055022020 1 



tion likelihood to 5. That yielded the detection of 130 
point sources as well as two regions of diffuse X-ray 
emission. 

A combined EPIC MOS1, MOS2, and PN im- 
age of Berkeley 87 and its surroundings is shown in 
Fig. [3] Two prominent hard (blue color) extended X- 
ray sources are located in the ON 2 SFR, slightly off 
to the north from the image center. The combined X- 
ray and IR images of ON 2 are shown in Fig. [4] The 
extended X-ray emission traces the eastern edge of 
ON 2S region and is also observed from the H n region 
G75.84+0.40inON2N. 

The angular resolution of XMM-Newton is 56". To 
exclude the potential confusion of extended emission 
with an unresolved stellar point source, we inspected 
optical and infrared images with higher angular reso- 
lution. The region around the Berkeley 87 cluster has 
been partially imaged in the mid-IR with the Sp itzer 
InfraRed Array Camera (IRAC: lFazio et al.l20o4 with 
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Fig. 4.— Image of ON 2 combined from XMM-Newton EPIC (0.25 to 12.0 keV band, blue) and Spitzer IRAC (3.6 fim, 
green, and 8.0//m, red). 



angular resolution of <2". The northern half of Berke- 
ley 87, including the H n regions GAL 75.84+0.40 and 
GAL 75.78+0.34, has observations in all four IRAC 
bands (3.6, 4.5, 5.8 and 8.0 fim) but the southern half 
of Berkeley 87 is only covered in 3.6 and 5.8 pm im- 
ages. We downloaded all applicable IRAC observa- 
tions from the Spitzer archive and combined the basic 
calibrated data to form mosaic images of the entire re- 
gion using the MOPEX software package. More in- 
formation on the instruments and pipeline processing 
can be found at the Spitzer Science Center's Observer 
Support website0. 



http://ssc. spitzer.caltech.edu/ost. 



5. Brief analysis of distribution of IR sources in 
Berkeley 87. 

YSOs can be identified by their IR excess, since 
they are still surrounded by dusty disks and envelopes 
that absorb stellar light and radiate at IR wavelengths. 
Theoretical predictions regarding the location of YSOs 
in color-magnitude diagrams (CMDs) and color-color 
diagrams (CCDs) can be used to compare with obser- 
vations in order to study the YSOs. Based on Spitzer 
IRAC and MI PS observations, this m ethod was re- 
cently used by iGruendl & Chu (l2009h to search for 
YSOs in the Large Magellanic Cloud. We use a sim- 
ilar method to search for objects with near- and mid- 
infrared excess in Berkeley 87. 
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Fig. 5. — CMDs and CCD based on IRAC and 2MASS photometry. The panels marked "all" use the photometric 
data for the entire field, while the panels marked "south" use a subset which exclude the ON 2 SFR, and the panels 
marked "north" correspond to the ON 2 SFR. In the leftmost panels of each row the arrow indicates the reddenning 
vector expected for Ay= 10 dCardelli et alJI 19891) and the solid lines in in the K s versus / - K s CMD and J — H versus 
H - K s CCD correspo nd to the expected location of the zero-age main sequence for an assumed distance of 1.23 kpc 
dBesseU & Brettlll988l) . Objects plotted with a filled red circle have a possible X-ray point source counterpart. Objects 
plotted with green circles are more than 1-c r beyond the black dashed line (/ - H = 1.692// - K) indicating possible 
excess infrared emission ( Allen et al1l2004 ). 




5 - 



q 
06 



10 - 




2 4 

[4.5]-[8.0] 



Fig. 6.— An [8.0] versus [4.5]-[8.0] CMD showing 
all sources in the Berkeley 87 region. Red points indi- 
cate sources from Figure|5]with [3.6]-[5.8]>1.0. Back- 
ground galaxies and AGN should genera lly fall below 



the d ashed line, [8.0] = 14-([4.5]-[8.0] dHarvev et al 
20061) . 



Fig. 7. — Combined Spitzer IRAC (3.6 //m green, 
4.5 jim red) and XMM-Newton EPIC (blue) image of 
Berkeley 87. The red regions mark the positions of 
IRAC sources with [3.6]-[5.8]>1.0 which may be 
YSOs. Image size is «22'x20'. North is up, and east 
is left. 



Aperture photometry was performed on the 3.6, 
4.5, 5.8, and 8.0 fim IRAC images to obtain mid-IR 
flux densities for sources throughout the Berkeley 87 
region. The results at each wavelength were com- 
bined for sources with <1'.'5 positional coincidence. 
We also obtained near-IR (JHK S ) flux densities from 
the Two Micron All Sky Survey Point Source Catalog 
(2MASS PSC; ISkrutskie etafl 120061) and combined 
these with the mid-IR measurements, again requiring 
a < 175 positional coincidence for a positive match be- 
tween sources detected in different bands. 

The resulting catalog does not have complete spa- 
tial coverage at all wavelengths over the entire re- 
gion discussed in this paper. Specifically there are no 
IRAC observations at 4.5 and 8.0 fim for declinations 
south of +37°20' (the southern and unobscured half of 
Berkeley 87). Therefore to search for evidence for a 
population of YSOs and PMSs we have confined our- 
selves to using the 2MASS JHK S bands and the IRAC 
3.6 and 5.8 fim bands. In Fig.[5]we present color-color 
and color-magnitude diagrams (CCDs and CMDs) that 
illustrate the difference in the stellar population we see 
when we compare: (1) the entire Berkeley 87 region, 
(2) the southern, unobscured, portion of Berkeley 87 
and (3) the northern portion of Berkeley 87 which 



contains the SFR ON 2. Specifically the sources for 
the entire Berkeley 87 region are within a box cen- 
tered at RA=20 h 21 m 37?7, DEC=+37° 24' 43'.'5, with 
ARA=20', ADEC=20', while the southern region 
is defined by a box centered at RA=20 h 21 m 37f7, 
DEC=+37° 17' 2275, with ARA=20', ADEC=10', (i.e. 
excluding the ON 2 SFR), and the northern region 
is defined by a box centered at RA=20 h 21 m 37?7, 
DEC=+37° 28' 2874, with ARA=20', ADEC=12', (i.e. 
the ON 2 SFR and surroundings). While the near-IR 
CCDs and CMDs constructed from 2MASS data ap- 
pear similar in the north and south, the [3.6] versus 
[3.6]-[5.8] CMD, constructed from the IRAC observa- 
tions, reveal a higher density of sources with red mid- 
IR colors (excess mid-IR emission). After correcting 
for the area we find that for [3.6]-[5.8]>1.0 there are 
~A times as many red sources in the northern region 
than in the south (see lower panel in Pig.[5j. 

The lack of MIPS 24 fim observations over much 
of the region, the incomplete IRAC coverage, and the 
limited sensitivity of 2MASS prevent a more detailed 
assessment as to the true nature of these red sources. 
Nevertheless, we speculate that these are YSOs associ- 
ated with Berkeley 87. For instance, if we select from 
the J - H versus H - K, CCD those 2MASS sources 
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Fig. 8.— Combined XMM-Newton MOS2 (0.25 to 12.0 keV band, blue) and Spitzer IRAC (3.6yum, green, and 8.0,um, 
red) image of ON 2S. The image size is as 6.5' x 6.5'. 



with near-IR execss, we find 49 sources (green circles 
in the middle panel of FigO. Only three of these ob- 
jects with near-IR excess are among the sources with 
mid-IR excess identified in the [3.6] versus [3.6]-[5.8] 
CMD. Thus most of the sources with mid-IR excess 
escape detection with 2MASS. 

We can, however, rule out that the sources with 
[3.6]-[5.8]>1.0 are dominated by background galax- 
ies. In Fig.|6]we present a [8.0] vs. [4.5]-[8.0] CMD 
for all sources with available data and find that most 
of them lie above the cutoff that has been used to sepa- 
rate bac kground galax i es and AGN in other population 



2006; Kozlowski & Kochanek 



studie s (iHarvey et al. 
2009h . Hence we believe that there are strong indica 



tions that the IRAC sources in Berkeley 87 with mid- 
IR excess are dominated by YSOs. In order to confirm 
and better quantify the YSO population throughout 



this region, deeper near-IR and more complete mid-IR 
observations are required. 

To better illustrate the distribution of these candi- 
date YSOs in ON 2 we plot the locations of all sources 
with [3.6]-[5.8]>1.0 in Fig. [7] We have also searched 
for possible X-ray counterparts to the 2MASS and 
IRAC points sources by comparing their positions and 
requiring a coincidence better than 3" for a positive 
match. A total of 47 matches were found but only one 
of those matches had a significant mid-IR excess (see 
Fig-El. The near-IR fluxes and colors of nearly all the 
sources with possible X-ray counterparts are consis- 
tent with normal main sequence stars. 

The very low X-ray detection rate of young stars 
in Berkeley 87 is not surprising. The limiting sen- 
sitivity of our XMM-Newton observations is Fx ~ 
1 x 10~ I4 erg s cirT 2 . Taking into account the aver- 
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Table 1 : The OB type stars detected by XMM-Newton and their parameters 



Star 


RA (J2000) 


DEC (J2000) 


Sp.type 




A^ H [10 21 cm- 2 ] 




HDE 229059 


20 h 21 m 


15*37 


37° 24' 31.3" 


Blla 


5.6 


8.6 


-7.5 


Berkeley 87-4 


20 h 21 m 


19*25 


37° 23' 24.3" 


B0.2III 


4.7 


7.4 


-7.1 


BD +36°4032 


20 h 21 m 


38 s .67 


37° 25' 15.5" 


08.5111 


5.1 


7.8 


-7.3 



cm 

2X10 31 — "- 1 



2 . At the 1.23 kpc, this cor- 
erg s _1 . Studies of low-mass 



age absorbing column in the direction of Berkeley 87, 
Nh * 8 x 10 21 cm~ 2 , the unabsorbed flux for a ther- 
mal source of X-ray emission with kTx = 0.8 keV is 
F x ^xlO^ergs- 1 
responds to Lx 
stars in the Orion Nebular Cluster shown that there is 
a corr elation between stellar age, mass, and X-ray ac- 
tivity dFlaccomio et al.ll2003l) . The X-ray luminosities 
of a few Myr old low- and solar-mass PMSs are lower 
than the detection limit of our observations. 

The CMDs shown in Fig.|5]provide one more indi- 
rect evidence that the SFR ON 2 and star cluster Berke- 
ley 87 are located at the same distance. The magni- 
tudes of sources located in the ON 2 region do not 
differ significantly from the magnitudes of sources in 
southern part of Berkeley 87. If ON 2 were 4 times 
more distant than Berkeley 87 (see discussion in Sec- 
tion ^, lb . one would expect many more faint sources 
in the right panels in Fig. [5] The spatial distribution of 
YSOs across t he entire ON 2N (F ig. [7]) confi rms earlier 



sugges tions of iTurner & Forb es ( 1982) and lDent et al 



dl988l) that star formation occurred nearly simultane- 
ously over this whole SFR. 

6. X-ray emission from massive stars in Berke- 
ley 87. 

In general, all massive stars with s pectral types ear- 
lier than B1.5 are X-ray emitters. ITurner & Forbes 



d 19821) identified 17 stars in Berkeley 87 as having 
early OB spectral types. Among them only three were 
unambiguously detected in our XMM-Newton observa- 
tion. This is because the large interstellar extinction in 
the direction of Berkeley 87 hampers the detection of 
soft X-ray sources such as OB stars. 

A combined X-ray and IR images of ON 2 and 
ON2S are shown in Figs.|4]and[8] The O-type giant 
BD +36°4032 (Berkeley 87-25) is located at the south- 
ern tip of the ON2S. Neighboring it are two B-type 
stars, separated by 8" and 1576 from BD +36°4032. 
Among these only the O star is detected. 

All spectra in this paper were analyzed using the 



xspec software ( Arnaudl 19961) . 

The X-ray spectrum of BD +36°4032 can be fitted 
with a thermal plasma model (apec) with temperature 
kTx = 0.6 + 0.1 keV, and a neutral hydrogen absorp- 
tion column A^h = (8 + 1) x 10 21 crrT 2 . The unabsorbed 
model flux is F x ss 2 x 10~ 13 ergcirT 2 s _1 , correspond- 
ing to Lx ~ 4 x 10 31 erg s _1 . This is the X-ray brightest 
massive star in Berkeley 87. 

Two B-type stars, HDE 229059 (Berkeley 87-3) and 
Berkeley 87-4, are also detected in our observation. 
Their X-ray spectra are fitted using a thermal plasma 
model with kTx = 0.6 + 0.2 keV corrected for inter- 
stellar column A^ H = (8 ± 1) x 10 21 crrT 2 . 

We use th e UBV photomet r y and spectral types 
determined in ITurner & Forbes! (119821) to derive stel- 
lar bolometric luminosities and color excesses Eb-v- 
The neutral hydrogen column densities are then es- 
timated using A^h _ 5. x 10 21 £b-v H-atoms cm -2 
jBohlin&Savagelll98lh and listed in TableHJfor the 
stars that are detected by XMM-Newton. The absorp- 
tion columns inferred from the analysis of the X-ray 
spectra are in good agreement with those from UBV 
photometry. 

The ratio of X-ray and bolometric luminosities for 
the detected OB stars (see Table[T|) appear to be slightly 
lower than the t ypical value of 10~ 7 for OB stars (e.g. 
Oskinoval 120051) . This is because the soft X-rays are 
missing due to the large interstellar absorption. Our 
crude one-temperature spectral models fitted to the ob- 
served spectra underestimate the contribution from the 
soft spectral range. Bearing this in mind, the level of 
X-ray emission from the OB stars in Berkeley 87 ap- 
pears to be normal for stars of these spectral types. 

7. X-ray emission from the Hn region 
GAL 75.84+0.40 in ON 2N 

In this section we will address the Hn regions 
in the northern part of ON 2 (see Fig.Q]). While 
no X-ray emission is detected from GAL 75.84+0.36, 
GAL 75.84+0.40 is a spectacular source of diffuse X- 
rays. We consider it in detail below. 
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Fig. 9. — Zoomed ON2N region from combined 
Spitzer IRAC 3.6, 4.5, and 8.0 /im image are shown 
as red, green and blue, respectively. Crosses mark the 
positions of source A (to the west) and source B (to the 
east). The image size is ~L5xl.'2 with north up and 
east to the left. 

7.1. Morphology of H n region GAL 75.84+0.40 



Matthews et al.l (119731) conducted a radio survey of 
ON 2. They noticed that GAL 75.84+0.40 has a com- 
plex morphology, and suggested that it consists of two 
compact Hn regions - the first is GAL 75. 84+0.40 A 
to the east, and the second is GAL 75. 84+0.40 B to 
the west. Within each compact H n region, there ex- 
ists a corresponding point source in the 2MASS PSC. 
These two IR point sources are separated by «;14'.'3 (or 
0.09 pc at d=1.23kpc). Their coordinates and K$ mag- 
nitudes are given in Table|2] Corresponding sources at 
these positions can be also seen in the Spitzer images 
(see Fig.|9]l. 

In the IRAC bands source A has flux densities of 
52.1+3.7 mJy and 44.1+5.0 mJy at 3.6 and 4.5 ^m, 
respectively ([3.6]^9.32mag and [4.5]-9.02mag), 
which are also consistent with an embedded 05-6V 
type star. On the other hand, for source B, we found 
a flux density of 194.2+7.3 mJy at 4.5 pm ([4.5]-7.4) 



Table 2: 2MASS point sources in GAL 75.84+0.40 





RA (J2000) 


DEC (J2000) 




A 


20 h 21 m 37 s :98 


+37° 31' 15.23" 


9.628 


B 


20 h 21 m 39 s .07 


+37° 31' 09.27" 


9.665 



Fig. 10.— Combined XMM-Newton MOS2 (blue) and 
Spitzer IRAC 3.6 pm (green) and 8.0 pm (red) images 
of ON 2N. The image size is ~ 5'5 x 5'5 with north up 
and east to the left. 



which suggests a mid-IR excess if the 2MASS /iT-band 
source corresponds to a early type star. Inspection of 
the 2MASS and IRAC images for sources A and B re- 
veal that both are amid complex diffuse emission but 
that while source A appears to be a single point source, 
source B appears elongated at all bands. We suggest 
that source B may be either a multiple or its flux mea- 
surements may suffer significant contamination from 
the surrounding diffuse emission. 

presented an alternative to 
view on the morphology of 



Garav et al 



Matthews et al 




GAL 75.84+0.40. Based on VLA radio maps they was 
suggested that GAL 75.84+0.40 is an inhomogeneous 
shell of gas with outer radius ~ 17", ionized by a sin- 
gle O-type central star. The highest resolution VLA 
2 cm radio continuum maps show the two component 
structure which can correspond to the bright rims of 
the shell. The "two" co mpact H n regions suggested by 
Matthews et al.l ( 119731) can be the brightest parts of this 
shell. Similar shell morphology, with brighten rims 
has been seen in planetary nebulae, post-AGB stars, 
and supernova remnants. 



11 





Fig. 11. — The same as in Fig.|9]but with contours of 
X-ray emission (contours are spaced logarithmically). 

7.2. X-ray emission from GAL 75.84+0.40 

A combined X-ray and IR image of the north- 
ern part of ON 2 is shown in Fig.[l0l while Fig.fTTI 
displays the IR Spitzer image overlaid with con- 
tours of the X-ray emission. The latter is ex- 
tended and fills a nearly circular region centered at 
20 h 21 m 38, +37°31'14"(J2000), slightly offset from 
the 2MASS source "B" (Table[2]i. The source detec- 
tion task indicates that the X-rays are diffuse. 

The extracted spectrum of the diffuse X-ray emis- 
sion in GAL 75.84+0.40 is shown in Fig.[[2] The Cash 
statistics was used to fit the spectrum by both thermal 
and by power-law models. The spectrum has a rela- 
tively low signal-to-noise ratio and can be equally well 
fitted either by thermal or by non-thermal emission. 
The best fit parameters are included in Table[4] The 
spectrum is heavily absorbed. The absorbing column 
is A^h ~ 4.5 x 10 22 cm" 2 , corresponding to Eb-v ~ 9. 

Turner & Forbes! (119821) determine from opti- 



cal the interstellar absorption in the direction of 
GAL 75.84+0.40 as E B -v=\H. Using N H = 5.0 x 
IO^-Eb-v cm" 2 , the column density is A^h = 8.5 x 
10 21 cm" 2 , i.e. a factor of five smaller than inferred 
from the X-ray spectral fits. Thus, we observe X-ray 
emission from one or more deeply embedded sources. 



200 - 



100 




o 



0.5 

log (Energy) [keV] 

Fig. 12. — Upper panel: Same as in Fig.[l0] The white 
circle shows the spectrum extraction region. Lower 
panel: The PN (red) and MOS2 (black) spectra of 
GAL 75.84+0.40. The best fit thermal plasma models 
plotted as solid lines for the corresponding detectors. 
The parameters of the model are given in Table[4] 



7.3. Discussion on the Origin of the X-ray Emis- 
sion from GAL 75.84+0.40 

The ionizing stars of Hn regions are intrinsic 
sources of X-ray emission. Assuming that stars 
GAL 75.84+0.40 contribute to the observed X-ray 
emission, we can obtain constraints on their stel- 
lar type. From the X-ray spectroscopy we deter- 
mined Eb-v ~ 9. Adopting Ry = 3.1, one ob- 
tains Ay - RyEs-v ~ 28 mag. The ratio be- 
tween visual and /T-band extinction is Ay/A K = 8.9 



dMoneti et al.ll2001l) . From the 2MASS Catalog, both 
stars in GAL 75.84+0.40 have w 9.6 mag, corre- 
sponding to Mk = - 4. Comparing with PoWR s tellar 
atmosphere models dHamann & Grafenen 120041) . this 
/iT-band absolute magnitude corresponds to roughly an 
05 type star. 

A number of earlier estimates for the ioniz- 
ing source of G AL 75 .84+0.40 exist. For exam- 
ple, iDent et al.l 01988b used the far-infrared flux 
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Table 3: Decomposing extended region filled with X-ray emission in ON 2S 



Region (as in FigO RA (J2000) DEC (J2000) Emission Origin 



OB 20 h 21 ra 38 +37° 25' 15" thermal stellar wind OB stars 

H 20 h 21 m 40 +37°25'34" ? star?UCHn? 

CN 20 h 21 m 43 +37°26'33" ? AGN? UCH n ?, star? 

R 20 h 21 m 41 +37°26'07" ? star?UCHn? 

D 20 h 21 m 41 +37° 26' 35" non-thermal? magnetic field? YSOs 



coordinates of the centers of the regions 



to derive the spectral type of the ionizing star in 
GAL 75.84+0.40 to be later th at Q9.7 assumin g a dis 



tance of 1 kpc. In contrast, iMatthews et al.l (119731) 



found that two 08V stars were needed to account 
for the radio continuum measurements of the H n re- 
gion after adoptin g a distance of 5.5 kpc. Finally, 
Garay et al.l (119931) determined that a single 06V star 
could explain their radio continuum measurements af- 
ter assuming a distance of 4.1 kpc. The only way to 
reconcile these earlier results with the 05 spectral type 
derived from the 2MASS magnitudes and X-ray spec- 
troscopy is to assume that the X-ray source is embed- 
ded deeper than the ionizing stars. Note that this con- 
clusion is independent of the adopted distance. 

Thus, there are good arguments that the ionizing 
stars in GAL 75.84+0.40 are not the main source of the 
observed X-ray emission: the X-ray source is deeper 
embedded, extended, and hard or non-thermal. 

The X-rays in GAL 75.84+0.40 fill a nearly circular 
area, with no obvious IR source being correlated with 
its center (see FigfTTli. One of the possibilities to ex- 
plain this can be an embedded cluster of young stars 
behind the H n region, which we observe in projection, 
and which is the source of X-ray emission. However, 
given the small linear extension of the X-ray source 
(0. 1 pc) and the character of its spectrum, a cluster of 
pre-main sequence stars seems to be an unlikely expla- 
nation. Albeit some contribution from point sources is 
possible, we believe that a truely diffuse component is 
present in GAL 75.84+0.40. 

Only a handful of star-forming regions with similar 
properties of X-ray emission is presently known (see 
Table|5]l. Besides the neighboring region ON 2S which 
we discuss in the next section, hard diffuse emission 
on similar spat ial scales is obse rved in the massive 
SFR RCW38 dWolketalJl2002b . An old shell-type 
SNR was considered as a possible explanation. How- 
ever, this explanation seems implausible in the case 
of GAL 75.84+0.40. Similar to the conclusions by 




Fig. 13. — The same as in Fig. [8] but with over-plotted 
contours (log scale) of X-ray emission. The red circles 
are around young stars identified from Spitzer photom- 
etry. The regions discussed in the text are identified by 
letters and arrows. The white bar in the low left corner 
has 1' length. 



Wolk et al.l d2002l) we suggest that invoking magnetic 
fields in the SFR may be required to understand their 
hard X-ray emission. We will return to this point when 
discussing the X-rays from ON 2S in Sect.[8] 

8. The ON 2S region 

In this and the following section we consider the 
southern ON2S region (Fig.Q]) and the X-ray emis- 
sion detected in its vicinity. Fig. [8] shows combined 
X-ray and IR images of ON 2S. In striking contrast to 
ON 2N, the extended X-ray emission in ON 2S does 
not coincide spatially with the UCHu regions and 
maser sources. Instead, the extended X-ray emission 
is observed to the east of the photo-dissociation region 
traced by the Spitzer images. 

The ON2S region is complex and compact, and 



13 



our XMM-Newton data are unable to fully resolve the 
individual components of the extended region. Nev- 
ertheless, the J4.5" angular resolution of the MOS2 
camera allows an investigation of the morphology of 
the hot gas. FigureQj] displays the combined Spitzer 
and XMM-Newton image overlaid with contours of X- 
ray emission. We distinguish five X-ray sources in 
ON 2S, corresponding to the local maxima in the con- 
tours. The coordinates of these sources are given in 
Table[3] The regions designated "CN", "R", "H", and 
"OB" were detected as point sources by the source 
detection software. The region "D" is detected as a 
diffuse source. In the following, we first analyze the 
point-sources, and then the diffuse X-ray emission. 

8.1. X-ray emission from point sources in ON 2S. 

Faint X-ray emission is detected south from 
Cygnus 2N (region marked CN in Fig.|T3Tl. Al- 
beit this emission appears extended in the adaptively 
smoothed image, the source detection software finds 
a point source at this location, coinciding within 
077 with the faint star USNO-B1.0 1274-0505781 
(Bl = 18.08 mag). No point source is detected at this 
position in any of the Spitzer IRAC channels. Simi- 
larly, no continuum radio sour ce at this position was 
found bv lShepherd et al.l dl997l) . 

The X-ray emission from the point source in CN 
is heavily absorbed, Nn>2 x 10 22 cirT 2 , and hard (see 
spectrum in panel 3 in Fig.fPfl). The absorbed flux is 
F x = (3.0+ 1.3) x 10 _14 ergcm" 2 s _1 . Because of its 
high absorption we exclude that the source is a fore- 
ground star. It could be either an embedded young star, 
or a background object (a Galactic star or an AGN). 
We consider the latter possibility as more likely, due to 
the lack of Spitzer or radio counterparts. However, on 
the basis of XMM-Newton observations we cannot rule 
out the interesting possibility that the X-ray emission 
south of Cygnus 2N originates from diffuse gas, while 
the close coincidence with a faint star is accidental. 

The region "R" is «8" away from the radio source 
GRS 75.77+0.34, and «6" away from a young star we 
identified from Spitzer images. It is plausible, that the 
radio source and the young star are, in fact, one object. 
It is also plausible to assume that the observed X-rays, 
at least partly, are associated with the star. It is inter- 
esting to note that a very red object seen only at 8/im 
in Spitzer images is located in the area R. 

The X-ray source in the region "H" (see Fig. [13]) is 
located some 10" away from the Maser 075.76+00.34 
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Fig. 14. — Upper panel: Same as in Fig. [8] The 
spectrum extraction regions are shown. Lower panel: 
XMM-Newton PN spectra of three regions in ON2S. 
The spectrum of the are containing OB -type stars is 
shown in plot 1 . The spectrum of whole ON 2S is 
shown in plot 2. The spectrum of region close to 
Cyg 2N is shown in plot 3. 



dHofner & Church welll[l996T) . and 6" away from the 
optical star Berkeley 87-83. The X-ray emission from 
H is too faint to constrain its morphology, and to ex- 
tract useful spectral information. Hence the nature of 
the X-ray emission from this area remains unclear. It 
could be due to unresolved young star(s) or a diffuse 
source. 

The X-ray emission from the region "OB" is dom- 
inated by O star BD +36°4032, and was discussed in 
Sect.|6] 

In Fig. [14] we compare the spectra extracted from 
the whole ON 2S region and from the smaller areas 
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Fig. 15. — Upper panel: Same as in Fig. [8] The spec- 
trum extraction region is shown. Lower panel: The 
background spectra are shown in the upper plot. XMM- 
Newton PN, MOS 1 and 2 spectra of diffuse emission 
in the region D are shown in the lower plot. The MOS2 
best fit model is shown as solid line, parameters of the 
models are given in Table[4] 



"CN" and "OB". It is apparent that the spectrum of 
the OB region is much less affected by absorption. 
The spectrum extracted from the large area is strongly 
absorbed and quite hard, showing that hard X-ray 
emission from embedded sources dominates the whole 
ON 2S region. 

8.2. Extended hard X-ray emission from ON 2S. 

In the XMM-Newton images the X-ray emission 
from the region "D" in Fig. [15] appears to be truely 
diffuse. We cannot fully exclude the presence of dis- 
crete sources embedded in the diffuse emission, al- 
beit we did not find any such sources in the avail- 
able multiwavelength catalogs, radio, IR and X-rays 
images. From analysis of the XMM-Newton images, 
we estimate that the diffuse X-ray emission occupies 



~3200arcsec 2 . At a distance of 1.23 kpc this corre- 
sponds to «;0.16pc 2 (2.25 pc 2 at a distance of 5.5 kpc 
respectively). 

The spectrum of diffuse emission is shown in 
Fig. [15] The inferred absorption column density 
strongly exceeds the one found in the direction of 
neighboring OB stars. The quality of the spectra is 
inadequate to rule out or confirm the presence of spec- 
tral lines. There is no indication of stron g Fe lines at 
6.4 ke V, as sometimes observed in SFRs dTakagi et al 
2002b. 



The spectrum is quite hard and can be fitted ei- 
ther with a thermal or an absorbed power-law model. 
The thermal plasma model fits best with a tempera- 
ture of 200 MK (Tabled. Such a high temperature 
of interstellar gas within a small scale SFR cannot be 
plausibly explained. An absorbed power-law model, 
K(E[keV]y y , where K is in photons keV _1 cm _2 s _1 at 
1 keV, fits best with y = 1.4 + 0.3. This is a plau- 
sible value for the synchrotron emission spectrum. It 
prompts us to suggest that the relatively small volume 
at the west side of the SFR ON 2S is may be filled by 
synchrotron emission. 

9. Discussion on the origin of the extended X-ray 
emission from ON2 S 

The diffuse thermal X-ray emission in SFRs may 
result from a collection of unresolved point sources, 
cluster winds, and wind-blown bubbles. Non-thermal 
emission requires the acceleration of electrons, for 
which the diffuse shock acceleration (DSA) mecha- 
nism is often invoked. In this section, we discuss the 
applicability of these mechanisms to the diffuse X-ray 
emission from ON 2S. 

9.1. Population of Point Sources 

Unresolved X-ray emission from many low-mass 
YSOs present in a massive SFR can mimic diffuse 
emission. An underlying population of low-mass stars 
was invoked to explain the hard extended X-ray emis- 
sion observed in the Orion and the O mega Nebula 
dTownslev et alj|2003l iGiidel et al.ll2008l). The known 



massive YSOs in ON 2 are 10000 yr old. At this age, 
low-mass YSOs are X-ray active, but may be unde- 
tected in 2MASS or IRAC images. From Fig.[5]we es- 
timate that we can see ZAMS down to early G stars for 
Ay = 0, but with just Ay = 10 this is already restricted 
to late F types. To uncover a population of low-mass 
YSOs in ON 2 high-quality, sensitive near-IR (JHK) 
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Table 4: Parameters of thermal and non-thermal plasma models used to fit the spectra of diffuse emission in 
GAL 75.84+0.40 and in ON2S (see FigsjjjHHj^ 



GAL 75.84+0.40 








Thermal plasma 


Non-thermal plasma 








tbabs*apec 


tbabs *powerlaw 


Nn 


[cm 




(4.0 ± 0.3) x 10 22 


(4.8 ± 0.6) x 10 22 








kT = 2.8 ±0.3keV 


y = 2.7+0.2 


F a 


[erg 


cm^ 2 s _I ] 


2.5 x 10" 13 




[erg 




6 x 10 31 


ON2S 








Thermal plasma 


Non-thermal plasma 








tbabs*apec 


tbabs *powerlaw 


Nn 


[cm 




(1.3 ± 0.2) x 10 22 


(1.2 + 0.3) x 10 22 








kT = 17.7 + 9.6 keV 


y = 1.4 + 0.3 


F a 


[erg 


cm -2 s _1 ] 


2 x 10" 13 


if 


[erg 


s" 1 ] 


4 x 10 31 



a absorbed 
b unabsorbed 



imaging is needed. 

Using results of iFlaccomio et al.l (120031) obtained 
from the study of stars in the Orion Nebula, we es- 
timate that the X-ray luminosity of a 10000 year 
old YSO with a mass of O.5SM/M ;l is L x ~ 4 x 
10 30 ergs _I . Therefore, the presence of just ten such 
objects can easily explain the observed luminosity of 
the extended X-ray emission in ON 2S. 

There are, however, arguments against interpreting 
the emission from ON2S as being solely due to an 
unresolved population of YSOs. Firstly, as was dis- 
cussed in Sect.|5] the Spitzer data confirm earlier sug- 
gestions that star formation occurred nearly simultane- 
ously over this whole ON 2 SFR. The observed X-ray 
emission is so hard, that even if it were present in the 
areas with higher absorbing column we would still be 
able to detect it. Yet, the patch of hard diffuse X-ray 
emission is spatially confined. There is no reason to 
expect that low-mass YSOs are strongly clustered at a 
location that is away from the higher-mass YSOs. 

Secondly, the spectrum from an unresolved popula- 
tion of YSOs would be much softer than what we ob- 



serve in ON 2S. iTownsley et al.l (120031) used Chandra 
observations of the Omega and the Rosetta Nebulae to 
obtain a cumulative spectrum of YSOs. The absorbing 
columns to YSOs in these nebulae are similar to those 



we derive in ON 2S. lTownsley et al.l (120031) found that 
the cumulative spectrum of point sources can be well 
fitted with thermal plasma models with temperatures 



not exceeding 3.1 keV. This is similar to what we infer 
from fitting the X-ray spectrum of GAL 75.84+0.40, 
but in ON 2S the temperatures are much higher (see 
Table©. 

Therefore, while it is possible that unidentified 
YSOs are embedded in ON 2, it appears that diffuse 
emission may also be present there. 

9.2. Shocked Stellar Wind 

The hot cluster winds filling the volumes of dense 
mass ive star clusters are drive n by stellar winds and 
SNe dChevalier & Cleggl[l985l) . A cluster wind can be 
expected fro m a 4-6 Myr old massive cluster, such as 
Berkeley 87 dOskinoval2005l) . 

In essence, the X-ray luminosity and temperature 
of a wind from a cluster with the radius R c \ depends 
on the input of mechanical energy and mass. The mul- 
tiwavelength observations of Berkeley 87 do not show 
any evidence of a recent SNR. Therefore, to crudely 
estimate the expected luminosity and temperature of 
the cluster wind from Berkeley 87, we consider only 
mechanic energy and mass input produced by stellar 
winds. In Berkeley 87 the wind energy production is 
dominated by the WR star WR 142. Recently, the stel- 
lar parameters of WR 142 were updated, and a mass- 
loss rate of 7 x lO~ 6 M yr~', and wind velocity of 
5500 km s _1 (or 400 0km s" 1 allowing fo r stellar ro- 



tation) were derived (lOskinova et al 
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To roughly estimate the energy input from other 
massive stars, we assume that there are * 30 such stars 
in Berkeley 87 (see Section[3]l. The wind velocities of 
O and early B stars are in the range «1500-3000km 
s _1 , while wind velocities of later B stars are in the 



range ^300-1000 km s _1 (lLamers & Cassinellill 19991) . 
Mass-loss diagnostics based on comprehensive stellar 
wind models that account for macro-clumping were 
recently used to infer empirical mass -loss rates of a 
few xlO ^Moyr" 1 for typica l O stars ( lOskinova et al 



20071: ISundqvist et al.l 120091) . The mass-loss rates of 



later B-t ype stars are at lea st one order of magnitude 



smaller dSearle et al.l 120081) . The wind-energy input 



from the late B and A type stars can be neglected. 

There are eight ear ly-type OB giants and super- 
giants in Berkeley 87 ( IMassev et al.l 1200 lb . We as- 
sume that each of them has a terminal wind velocity of 
2000 km s _1 and a mass-loss rate of 10~ 6 M Q yr _1 . Fur- 
thermore, we assume that there are 15 stars with termi- 
nal wind velocities of 600 km s and mass-loss rates 
of lO _8 M yr _1 (corresponding to later B type stars). 
The cluster radius is 3 pc (Fig.|TJ. 

Usi ng scaling relations of IStevens & Hartwell 
(12003b . the cluster wind luminosity is Lx ~ 
10 M^(l? c iV*) ergs -1 , and the cluster wind temper- 

- 2 

ature is Tx ~ 15 V* MK (accounting for the mass- 
loading from "cool" matter can lead to the somewhat 
lower temperature.). The dimension of M is M yr _I , 
R c \ is in pc, and V» is in kms -1 . M* is the sum of stel- 
lar mass-loss rates in the cluster, and V* is the mean 
stellar wind velocity weighted by mass-loss rates. Us- 
ing the wind parameters of the stars in Berkeley 87, a 
very hot (T x ~ 100 MK) but faint (L x ;l0 30 erg s" 1 ) 
cluster wind emission can be expected. Note, that our 
estimates are not sensitive to the assumed wind param- 
eters of the OB stars, because the kinetic energy of the 
WO wind alone is an order of magnitude higher that 
the combined kinetic energy input from all other stars. 

Our XMM-Newton data do not show diffuse emis- 
sion filling the cluster, with temperature peaking in the 
cluster center, and extending beyond its boarders, as 
predicted by the analytical models. Instead, we ob- 
serve localized, small (;0.15pc 2 ) patches of diffuse 
emission within the cluster. Therefore, it is not possi- 
ble to attribute the diffuse X-ray emission from ON 2S 
to the cluster wind. 

If indeed the molecular cloud Onsala 2C and the ad- 
jacent star-forming region ON 2 are immediate neigh- 
bors of the cluster Berkeley 87, whose intra-cluster 



medium is filled with hot tenuous wind, the situ- 
ation may be analogou s to the one considered by 



Cowie & McKed (119771) for the evaporation of a cool 



cloud i n a hot gas. For a sphe rical cloud, the solu- 
tions of ICowie & McKed d 1977b predict that the tem- 
perature rises steadily with the distance r from the 
cold cloud (T oc (1 - r~ l ) 0A ) and the X-ray emis- 
sion forms a "halo" around t he cool cloud (see also 
Harper-Clark & Murrayf2009l) . Broadly speaking, the 



thermal conduction would lower the temperature and 
increase the X-ray luminosity o f the hot gas close to 
the interface ( St effen et al . 200§. 



There are, however, a number of arguments against 
attributing the X-ray emission from ON 2S to the ther- 
mal evaporation of cool cloud material. The observed 
X-ray temperature at the interface between cool and 
hot gas appears to be unreasonably high (~ 200 MK). 
Moreover, the diffuse X-ray emission is localized in 
small areas, while the interface between cool and hot 
gas in Berkeley 87 extends over parsecs (see sketch in 
Fig.©. 



Furthermore, as noticed e.g. by I Wang et al.l (12006b . 
the presence of the Fe-line complex at « 6.4 - 6.7 keV 
is a good indicator that cool clo ud material is in- 
volved in the generation of X-rays. I Wang et al. I d2006l) 
propose a cluster-cloud collision scenario, where the 
emission in Fe lines traces the shocked cloud gas. This 
line complex is very prominent in the Arches and the 
Quintuplet cluster spectra. Using the spe ctral param 



eters o f the Fe-line complex obtained in I Wang et al 



(12006b we conclude that if Fe emission of comparable 
strength were present in ON 2S it would be noticeable 
even in our low S/N spectra. However, there is no in- 
dication for the presence of the Fe line complex in the 
X-ray spectra of ON 2S. 

Thus, taking into account i) the spatial distribu- 
tion of diffuse emission, ii) the very high temperature 
needed to fit a thermal model to the observed spectrum, 
and Hi) the absence of fluorescent iron lines in the 
spectra, we conclude that the diffuse X-ray emission 
is not a consequence of the interaction of the cluster 
wind from Berkeley 87 with the cool cloud Onsala 2C. 

Three massive stars (two BI and one OIII) are lo- 
cated just south from the area filled with extended X- 
ray emission (Fig.[8j. The local interaction between 
the winds of these three stars may, in principle, lead 
to the heating of nearby regions. This scenario, how- 
ever, is not confirmed by the data. As discussed in Sec- 
tion l8.ll the X-ray luminosity and spectrum extracted 
from the region around these three OB stars is domi- 
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nated by X-ray emission of BD+36°4032 and does not 
require any additional sources such as interactions be- 
tween stellar winds. 



Polcaro et al 



( 11991b reported the detection of 
strong diffuse emission around WR 142 in the optical. 
They argue that this diffuse emission likely originates 
in a supersonic flow centered on WR 142, and can be 
considered as evidence of a hot bubble around this star. 
Although theoretically expected, there is a dearth of 
detected diffuse X-ray e mission from wind-blown bub- 



bles a round WR stars dChu et alj l2003t IWrigge et al 



2005). The only two detected hot bubbles show a limb- 



brightened morphology and are extended on the scale 
of parsecs. The size of a hot bubble around WR 142, 



estima ted using the classical work by IWeaver et al 
dl977l) . should be « 40 pc or Si. 5° at the distance 
of Berkeley 87. It should be noted, however, that 



IWeaver et al.l ( 1977 ) considered the case of an iso- 
lated star in a uniform medium, i.e. conditions that 
are clearly not valid for WR 142. A superbubble and 
super-sh ells can be expe cted around rich OB asso- 
ciations (IChu et al.l 120031) . such as Cyg OB 1, which 
comprises star clusters Berkeley 86, 87, IC 4996, and 
NGC6913. However, these structures are extended 
over large scales of ~ 10 2 pc. The location of Berke- 
ley 87 in the complex region of the Cygnus X super- 
bubble makes the detection of an individual bubble 
around Berkeley 87 even more difficult. 

Thus we conclude that stellar winds cannot be the 
main reason for the observed apparently diffuse emis- 
sion. 

9.3. Synchrotron Emission in Random Magnetic 
Fields 

Diffuse hard X-ray emission with a power-law spec- 
trum can be produced by the synchrotron mechanism. 
The latter requires the presence of energetic particles 
and magnetic fields. 

It has been argued that the shocks in the clus- 
ter winds and/or colliding wind binaries efficiently 
acceler ate electrons and p rotons to relativistic en- 
ergies ( Bvkov et al. 200 It 



20011: lOuataert & Loebl 12005 



Pittard & Dougherty 20061) . The model calculations 
by lBednarek (2007) show that particles up to TeV en- 
ergies should be present in Berkeley 87. These parti- 
cles should be advected from the cluster on the time 
scale of ~ 10 3 yr. A simple scaling relations to es- 
timate the expected surface brightness of the syn- 
chrotron emission from the cluster of massive stars 



were derived by Ouataert & Loebl ( 2005 ) and applied 
to the cluster in the central parsec of the Galaxy. 
Berkeley 87 is less rich than the Galactic Center clus- 
ter, but it contains WR 142, a WO star that alone has 
the wind kinetic energy « 1.5 x 10 38 erg s , making it 
comparable to the whole of the Galactic Center clus- 
ter. The small number of OB stars in Berkeley 87 re- 
sults in the order of two lower UV luminosity, thus 
the cooling of particles by inverse Compton process 
is less significant. Accord ing to the scaling relations 
of lOuataert & Loebl d2005l) . the surface brightness of 
Berkeley 87 in y-rays is ~ 100 lower than that of the 
Central Cluster. 

It is interesting to note that the winds from the dense 
group of three OB stars south of the ON 2S (see sketch 
in Fig. [2]) may provide an additional source of particles 
accelerated in situ nearby to ON 2S. 

To estimate at which distance from the stars the 
winds terminate and shock the Hn gas, we consider 
a typical O star wind density profile which declines 
with distance as r~ 2 . The wind density close to the sur- 
face of an O star is ~ 10 10 cirT 3 . For the typical ambi- 
ent matter density, ~ 1 cirT 3 , the stellar wind pressure 
would be equal to the ambient matter pressure roughly 
at the distance of the patch of diffuse X-ray emission 
(i.e 0.3 pc away from OIII star BD+36°4032). 



Recently, iBykov et al.1 (120081) addressed the effect 
of a random magnetic field on synchrotron emission. 
In was shown (in the context of young supernova rem- 
nants) that prominent localized structures can appear 
in synchrotron maps of extended sources with ran- 
dom magnetic fields, even if the particle distribution is 
smooth. The bright structures originate as high-energy 
electrons radiate efficiently in local enhancements of 
the magnetic field. The size of the "patch" of hard X- 
ray emission in ON 2S is about 0. 1 pc 2 . This could be 
relevant in the framework of the Bykov et al. model 
if the shock velocity is a fewxl0 3 kms~' (Bykov A., 
priv. communication). This is a plausible number for 
a stellar as well as for a cluster wind. 

Magnetic fields, along with turbulence, play a lead- 
ing role in star formation ( Crutcher et al. 20091) . The 
field strengths (50-700 pG) of an order, or even two, 
higher than the field in the diffuse ISM are detected 
in SFRs, moreove r the field geomet ry is shown to be 
far from uniform d S chleunin dl 1 99 8l) . These measure- 
ments provide strong support to the idea that magnetic 
fields of similar strengths can be present in other, if 
not all , SFRs, with ON 2S being no exception. iFerland 
d2009h reviewed several recent observational studies of 
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Table 5: High-mass SFRs where hard diffuse X-ray emission was detected 3 



1? m r\r\ 






Size 


Tl1 Ctfltl/^f 1 


Model 


V^UllllllCllL 


T? pfprptipp 

JXCICICIIA-C 




[10 32 erg s _I ] 


[10 22 cm~ 2 ] 


[pc] 


[kpc] 


thermal b 
[keV] 


power law 

r 






Orion 


0.6 


0.1 


2 


2.6 


0.4 




point sources (?) 


(1) 


ON2S 


0.4 


1.5 


0.1 


1.2 


18 


~ 1.4 


synchrotron (?), point sources 


(2) 


GAL 75.84+0.40 


0.6 


4.6 


0.1 


1.2 


3 


* 2.6 


point sources (?) 


(2) 


RCW38 


1 


1 


1.5 


x2 


1.7 


* 2.8 




(3) 


NGC6334 


0.1-5 


0.5-10 


I 


1.7 


> 1 






(4) 


Westerlund 1 


300 


2 


10 


5 


S3 




IC (?) 


(5) 


NGC 3603 


200 


0.7 


4 


7 


3 






(6) 


Arches 


200 


10 


3 


S.5 


5.7 






(7) 


Sgr B2 


9 


40 


0.2 


8.5 


10 






(8) 


W49A 


30 


50 


0.3 


11.4 


7 






(9) 


LMC 30 Dor C 


10 4 


0.1 


10 


50 




* 2.5 


superbubble 


(10) 


LMCN11 


1.5 x 10 3 


0.5 


10 


50 




a 1.7 


superbubble 


(11) 


LMC N51D 


4x 10 3 


0.03 




50 




as 1.3 


superbubble 


(12) 



"All numbers in this Table are approximate. The reader is urged to consult the original publications describing these complex 
objects and their observations in detail. 

b Temperature of the hardest component in multi-temperature spectral fits 
c Comment on the properties of suspected non-thermal X-ray emission, if any 

(1) Gudel et al. (2008) ; (2) this work; (3) Wolketal. (2002) : (4) Ezoe et al. (2006) ; (5) Muno (2003) ; 

(6) Moffat et al. C2002) : (7) Yusef-Zadeh et al. f2002) : (8) Takagi et al. f2002) ; t9) Tsuiimoto' et al. f2006) : 

(10) Bamba et al. (2004) ; (11) Maddox et al. (2009) ; (12) Cooper et al. (2004) ; 



the relationships between magnetic fields, stellar feed- 
back, and the geometry of Hn regions. The observa- 
tions reveal that magnetic field lines can be preferen- 
tially aligned perpendicular to the long axis of the qui- 
escent cloud before stars form. After star formation 
and push-back occurs, ionized gas will be constrained 
to flow along the field lines and escapes from the sys- 
tem in directions perpendicular to the long axis. Wave 
motions may be associated with the field and so could 
contribute a turbulent component to the observed line 
profiles. 

It appears that all ingredients required to produce 
a small bright patch of synchrotron radiation in the 
X-ray image may be present in ON2S. The particles 
accelerated up to TeV energies may result from the 
shocked cluster winds and/or from the YSO outflows. 
Sufficiently strong turbulent magnetic fields may be 
present in the vicinity of star forming regions. In addi- 
tion, the presence of clumps of matter w ould enhance 



the s urface brightness of the radiation dBvkov et al 



2008), and such clump s of matter are presen t in the ra- 
dio maps of the region ( Shepherd et al.fl997 ). To sum- 
marize, a synchrotron emission may be expected from 
the region where the strongly shocked stellar winds 
and turbulent magnetic fields co-exist. 



10. Comparison with other SFRs 

In Table[5] we exp and the list compiled by 



Tsuiimoto' et alj (12006b to summarize the properties 
of high-mass star forming regions with detected hard 
diffuse X-ray emission. 

It can be immediately seen that the properties of 
SFRs in Tableware quite diverse. It appears that three 
broad categories can be distinguished. The superbub- 
bles around large clusters and associations of mas- 
sive stars belong to one of these categories. Large 
spatial scale (10-100pc), high luminosity, and often 
limb brighten morphology are characteristic for ob- 
jects such as 30DorC, LMCN11, LMCN51D. The 
star clusters blowing cluster winds belong to the sec- 
ond category. The examples include Westerlund 1, 
the Arches, and NGC 3603. Although their X-ray 
emission is hard, it can, usually be fitted with hot 
(r x 530MK) thermal plasma. The diffuse X-rays fill 
the cluster interior on the scale of a few pc. In the third 
group, where, as we believe ON 2 belongs, small scale 
areas (~ 0.1 pc) are filled with diffuse X-rays. In such 
SFRs as ON 2, RCW38, WR 49 A, or Sgr 2, the hot 
(JxSlOOMK) or non-thermal plasma is found in im 



media te vic inity of UCH n regio ns. While iTakagi et al 



(12002b and iGudel et al.1 d2008l) suggest that hard X- 
rays may originate from unresolved populations of 
low-mass pre-main sequence stars, we suggest, that at 
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least in some objects, the non-thermal emission may 
be expected from interactions between magnetic fields 
and the particle accelerated in shocked stellar winds. 

It is interesting to note that the distribution of hard 
extended X-ray emission in the two closest SFRs - 
Orion nebula and ON 2S - is somewhat similar: diffuse 
X-ray emission filling a cavity seen in the IR images, 
slightly offset from a small group of massive stars. The 
X-ray spectrum, however is noticeably harder in the 
case of ON2S. 

11. Concluding remarks 

Throughout the paper we adopted the distance d = 
1 .23 kpc to both the young massive star cluster Berke- 
ley 87 and the SFR ON 2. This assumption has allowed 
us to propose physical interaction between the radia- 
tive and mechanical feedback from OB stars in Berke- 
ley 87 and the observed diffuse X-ray emission in the 
regions of active star formation. 

In the majority of papers published on the SFR 
ON 2, a distance in excess of 4 kpc is adopted. In this 
case, there would be no physical connection between 
the massive star cluster and the SFR. The X-ray lumi- 
nosities and the areas filled with X-ray emission would 
be factor of ;16 larger. The results of our spectral anal- 
ysis would, however, not been affected. 

To summarize, we have conducted XMM-Newton 
observations of the massive star-forming region ON 2. 
The observations and their subsequent analysis have 
shown: 

1. Diffuse X-ray emission on a cluster scale, which 
might have been expected to result from the cluster 
wind and wind blown bubble in Berkeley 87, was not 
detected. 

2. From a literature search we do not confirm previous 
reports of y-ray emission from Berkeley 87. 

3. The northern (ON2N) and the southern (ON2S) 
parts of the star forming complex ON 2 are bright 
sources of diffuse X-ray emission. 

4. X-ray images of ON 2N show that X-ray emis- 
sion fills the interior of the compact Hn region 
GAL 75.84+0.40. This emission is diffuse and 
strongly absorbed, indicating that it originates from 
deeply embedded sources. We rule out the ionizing 
stars of the H n region as its origin, and speculate that 
it can result from an embedded cluster of young low 
mass stars. 



5. In ON2S the extended X-ray emission traces the 
eastern edge of this SFR. This extended emission con- 
sists of point sources superimposed on or immersed in 
diffuse emission. 

6. We discuss different possible scenarios to explain 
the apparently diffuse emisison from ON 2S, such as 
cluster wind, synchrotron radiation, and unresolved 
point sources. We favor the last two options as the 
most probable explanations. 

Note added in proof. New Chandra X-ray telescope 
observations of a part of Berkeley 87 became public 
after this article had been accepted. The Chandra 
exposure is less sensitive than the XMM-Newton ob- 
servations discussed in the present paper, but has a 
superior angular resolution of il". The new Chan- 
dra images reveal three discrete sources in the region 
termed D in Figs. [13] [15] According to a first anal- 
ysis, the Chandra spectra of these three point sources 
correspond to strongly absorbed thermal emission with 
N H ~ 3 x 10 22 cirT 2 and kT « 2keV. Such param- 
eters are usual for YSOs as discussed in Section l9TI 
Hence the new data indicate that unresolved point 
sources contribute to the apparently diffuse emission 
from ON2S detected by XMM-Newton, confirming 
what we have discussed in Section[9] as one possibil- 
ity. However, the sum of the fluxes from the discrete 
sources is about three times smaller than the flux from 
the whole region D as measured with XMM-Newton, 
still leaving open the question about the true nature 
of the major part of extended emission observed with 
XMM-Newton. 
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A. On the absence of detected y-ray emission from Berkeley 87 



It has long been realized that young massive stellar associations can be potential sources of cosmic- and y-rays, 
where the acceleration of particles occurs in the shocks produced by t urbulent interactions of supersonic stellar winds 
dMontmerlelll979t ICasse & Pau3ll980t ICesarskv & Montmerielll983b . Berkeley 87 is a testbed for the models of y- 
ray production due to its proximity, relatively high space density of massive stars, and, chiefly, due to the cluster 
membership of WR 142 - the star with perhaps the most powerful stellar wind in the Galaxy. 

However, our careful check revealed that in contrast to the reports in the literature, the high-energy emission from 
Berkeley 87 was not previously observed. 

We were not abl e to confirm the detection of diffuse X-ray emission from Berkeley 87 by EXOSAT as referenced 
bv lBednarekl d2007l) . Berkeley 87 was observed ~ 1 deg off-axis. There is no reference in the literature to the detection 
of Berkeley 87 or any object within its boarders by EXOSAT. Similarly, no detection of Berkeley 87 by ASCA was 



ever reported. [Roberts et al.l d2002l) report the identification of an ASCA X-ray source coincident with the EGRET 



y-ray source 2CG 075+00, GeV J2020+3658 or 3EG J2021+3716 as the young radio pulsar PSR J202 1+3651, which 
is located at a distance of * lOkpc. The pulsar and its associated pulsar wind nebula G75. 2+0.1 are w 28' away from 
the outer border of Berkeley 87. 



Bednarekl (120071) considered two EGRET sources, 3EG J2016+3657 and "3EG J2021+4716" , as related to Berke- 



ley 87_£hj^ter L The_elliptical fits to the 95% confidence position contours of all 3EG EGRET sources are presented 
in iMattox et al.l (1200 lb . The smallest separation of the outer border of Berkeley 87 from 95% confi dence position 
contours of 3EG J2016+3657 is as 30'. No source "3EG J2021+4716" is present in the 3EG catalog dHartman et al. 



1999b . 



Aharonian et al. d2006b report that the observations of Berkeley 87 region yielded only uppe r limits that are a 
factor of <3 lower than TeV fluxes expected from an OB stellar association. The Milagro group dAbdo et alj|2007b 



reported the detection of an extended TeV source MGRO J2019+37, but their analysis favors its associations with the 
PSR J2021+3651 and its pulsar wind nebula G75.2+0.1. 

From our literature search, we conclude that prior to our XMM-Newton observations, there were no firm detection 
of X-ray diffuse sources from Berkeley 87. We also conclude that there are no y-ray source identifications within 
Berkeley 87. 



22 



REFERENCES 

Abdo, A. A. et al. 2007, ApJ, 658, L33 

Aharonian, F. et al. 2006, A&A, 454, 775 

Allen, L. E., et al. 2004, ApJS, 154, 363 

Arnaud, K. A. 1996, ASP Conf. Ser. Vol. 101. Astron. 
Soc. Pac, San Francisco, p. 17 

Bamba, A. et al. 2004, ApJ, 602, 257 

Bohlin, R. C. & Savage, B. D. 1981, ApJ, 249, 109 

Bednarek, W. 2007, MNRAS, 382, 367 

Bessell, M. S., & Brett, J. M. 1988, PASP, 100, 1134 

Brogan, C. L. & Troland, T. H. 2001, ApJ, 560, 821 

Bykov, A. 2001, Space Sci. Rev., 99, 317 

Bykov, A. M., Uvarov, Y. A., & Ellison, D. C. 2008, 
ApJ, 689,L133 

Cardelli, J.A., Clayton, G.C, & Mathis, J.S. 1989, 
ApJ, 345, 245 

Casse, M. & Paul, JA. 1980, ApJ, 237, 236 

Cesarsky, C.J. & Montmerle, T. 1983, Space Sci. Rev., 
36, 173 

Chevalier, R. A. & Clegg, A. 1985, Nat, 317, 44 

Chu, Y.-H., Guerrero, M. A., Gruendl, R. A., Garcia- 
Segura, G, & Wendker, H.L. 2003, ApJ, 599, 1 189 

Cooper, R. et al. 2004, ApJ, 605, 751 

Cowie, L. L. & McKee, C. F. 1977, ApJ, 21 1, 135 

Crutcher, R. M., Hakobian N, & Troland, T. H. 2009, 
ApJ, 692, 844 

Dent, W.R.F, Macdonald, G.H., & Andersson, M. 
1988, MNRAS 235, 1397 

Ezoe, Y, Kokubun, M., Makishima, K., Sekimoto, Y, 
& Matsuzaki, K. 2006, ApJ, 638, 860 

Fazio, G, et al. 2004, ApJS, 154, 10 

Ferland, G. J. 2009, IAUS, 259, 25 

Flaccomio, E., Damiani, F, Micela, G, Sciortino, S., 
Harnden, F. R. Jr, Murray, S. S., & Wolk, S. J. 2003, 
ApJ, 582, 398 



Garay, G, Rodriguez, L.F., Moran, J.M., & Church- 
well, E. 1993, ApJ, 418, 368 

Gruendl, R. A., & Chu, Y.-H. 2009, ApJS, 184, 172 

Gudel, M., Briggs, K.R., Montmerle, T., Audard, M., 
Rebull, L., & Skinner, S.L. 2008, Science, 319, 309 

Hamann, W.-R. & Griifener, G. 2004, A&A, 427, 697 

Harper-Clark, E. & Murray, N. 2009, ApJ, 693, 1696 

Hartman, R.C., et al. 1999, ApJS, 123, 79 

Harvey, P. M., et al. 2006, ApJ, 644, 307 

Hofner, P. & Churchwell, E. 1996, A&AS, 120, 283 

Kozlowski, S., & Kochanek, C. S. 2009, ApJ, 701, 508 

Maddox, LA., Williams, R.M., Dunne, B.C., & Chu, 
Y.-H 2009, ApJ, 699,911 

Massey, P., DeGioia-Eastwood, K., & Waterhouse, E. 
2001, AJ, 121, 1050 

Mathys, G. 1987, A&AS, 71, 201 

Matthews, H.E., Goss, W.M., Winnberg, A., & 
Habing, H.J. 1973, A&A, 29, 309 

Matthews, N, Anderson, M., & MacDonald, GH. 
1986, A&A, 155, 99 

Mattox, J., Hartman, R.C., & Reimer, O. 2001, ApJS, 
135, 155 

Meynet, G. & Maeder, A. 2005, A&A, 429, 581 

Moffat, A.F.J., et al. 2002, ApJ, 573, 191 

Moneti, A., Stolovy, S., Bommart, J. A. D. L, et al. 
2001, A&A, 366, 106 

Montmerle, T., 1979, ApJ, 231, 95 

Muno, M.P, et al. 2003, ApJ, 589, 225 

Negueruela, I. 2004, AN, 325, 380 

Lamers, H. J. G. L. M. & Cassinelli, J. P. 1999, Intro- 
duction to Stellar Winds (Cambridge: Cambridge 
Univ. Press) 

Levesque, E.M., Massey, P., Olsen, K.A.G, Plez, B., 
Josselin, E., Maeder, A., & Meynet, G, 2005, ApJ, 
628, 973 

Oskinova, L.M. 2005, MNRAS, 361, 679 



23 



Oskinova, L. M., Hamann, W. R., & Feldmeier, A. 

2007, A&A, 476, 1331 

Oskinova, L.M., Hamann, W.-R., Feldmeier, A., Ig- 
nace, R., & Chu, Y.-H. 2009, ApJ, 693, L44 

Palagi, F., Cesaroni, R., Comoretto, G., Felli, M., Na- 
tale, V., 1993 A&AS, 101, 153 

Pittard, J. M. & Dougherty, S. M. 2006, MNRAS, 372, 
801 

Polcaro, V.F, Giovannelli, F, Manchanda, R. K., 
Norci, L., & Rossi, C. 1991, A&A, 252, 590 

Pravdo, S.H., Tsuboi, Y., Suzuki, Y, Thompson, T.J., 
Rebull, L. 2009, ApJ, 690, 850 

Quataert, E. & Loeb, A., 2005, ApJ, 635, L45 

Reifenstein, E.C., Wilson, T.L., Burke, B.F, Mezger, 
P.G., & Altenhoff, W.J. 1970, A&A, 4, 357 

Roberts, M.S.E. et al. 2002, ApJ, 577, L19 

Sundqvist, J. O., Puis, J., & Feldmeier, A. 2009, arXiv 
0911.3361S 

Schaller G, Schaerer D., Meynet G, Maeder A. 1992, 
A&AS, 96, 269 

Schneider, N., Simon, R., Bontemps, S., Comeron, E, 
& Motte, F. 2007, A&A 474, 873 

Schleuning, D. A. 1998, ApJ, 493, 811 

Searle, S. C, Prinja, R. K., Massa, D., & Ryans, R. 

2008, A&A, 481,777 

Shepherd, D.S., Churchwell, E., & Wilner, D.J. 1997, 
ApJ, 482, 355 

Skrutskie, M. E, et al. 2006, AJ, 131, 1163 

Steffen, M., Schonberner, D., Warmuth, A. 2008, 
A&A, 489, 173 

Stevens, I.R. & Hartwell, J.M. 2003, MNRAS, 339, 
280 

Takagi, S., Murakami, H., & Koyama, K. 2002, ApJ, 
573, 275 

Townsley, L.K., Feigelson, E.D., Montmerle, T., 
Broos, PS., Chu, Y.-H., & Garmire, GP 2003, ApJ, 
593, 874 



Tsujimoto, M., Hosokawa, T., Feigelson, E.D., Get- 
man, K.V., & Broos, PS. 2006, ApJ, 653, 409 

Turner, D. G. & Forbes, D. 1982, PASP, 94, 789 

Turner, D. G, Rohanizadegan, M., Berdnikov, L. N., 
& Pastukhova, E.N. 2006, PASP, 118, 1533 

Uyaniker, B., Furst, E., Reich, W., Aschenbach, B., & 
Wielebinski R. 2001, A&A, 371, 675 

Vanbeveren, D., de Donder, E., van Bever, J., van 
Rensbergen, W., & de Loore C. 1998, New Astron., 
3,443 

Weaver, R., McCray, R., Castor, J., Shapiro, P., & 
Moore, R. 1977, ApJ, 218, 377 

Wang, Q. D., Dong, H, & Lang, C. 2006, MNRAS, 
371,38 

Wink, J. E., Altenhoff, W. J., & Mezger, PG 1982, 
A&A, 108, 227 

Wolk S. J., Bourke T. L., Smith R. K., Spitzbart B., & 
Alves, J., 2002, ApJ, 580, L161 

Wrigge, M., Chu, Y.-H., Magnier, E.A., & Wendker, 
H.J. 2005, ApJ, 633, 248 

Yusef-Zadeh, E, Law, C, Wardle, M., Wang, Q. D., 
Fruscione, A., Lang, C. C, & Cotera, A. 2002, ApJ, 
570, 665 



This 2-column preprint was prepared with the AAS IATgX macros 
V5.2. 



24 



